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ABSTRACT 

Context. It has been established that Gamma-Ray Bursts (GRBs) are connected to Supernovae (SNe) explosions of Type Ib/c. 
Aims. We intend to test whether the hypothesis of Type Ib/c SNe from different massive progenitors can reproduce the local GRB rate 
as well as the GRB rate as a function of redshift. We aim to predict the GRB rate at very high redshift under different assumptions 
about galaxy formation and star formation histories in galaxies. 

Methods. We assume different star formation histories in galaxies of different morphological type: ellipticals, spirals and irregulars, 
which have been already tested in self-consistent galaxy models reproducing both chemical and photometrical properties of galaxies. 
We explore different hypotheses concerning the progenitors of Type Ib/c SNe: i) single massive stars (M > 25M e , Wolf-Rayet stars), 
ii) massive close binaries (12-20M o ) and iii) both Wolf-Rayet stars and massive binaries. We conclude that the mixed scenario (iii) is 
preferable to reproduce the local Type Ib/c SN rates in galaxies and we adopt this scenario for comparison with the GRB rates. 
Results. We find an excellent agreement between the observed GRB local rate and the predicted Type Ib/c SN rate in irregular 
galaxies, when a range for single Wolf-Rayet stars of 40-100M o is adopted. We also predict the cosmic Type Ib/c SN rate by taking 
into account all the galaxy types in an unitary volume of the Universe and we compare it with the observed cosmic GRB rate as a 
function of redshift. By assuming the formation of spheroids at high redshift, we predict a cosmic Type Ib/c SN rate, which is always 
higher than the GRB rate, suggesting that only a small fraction (0.1-1%) of Type Ib/c SNe become GRBs. In particular, we find a 
ratio between the cosmic GRB rate and the cosmic Type Ib/c rate in the range 1CT 2 - 1(T 3 , in agreement with previous estimates. 
Finally, due to the high star formation in spheroids at high redshift, which is our preferred scenario for galaxy formation, we predict 
more GRBs at high redshift than in the hierarchical scenario for galaxy formation, a prediction which awaits to be proven by future 
observations. 

Key words. Galaxies: evolution; Galaxies: fundamental parameters; Supernovae: general; Gamma-Rays: bursts. 



. ^ ' 1. Introduction lines exhibiting P-Cygni profiles and represent about 70% of the 

k> , s , , , r , r, , r exploding stars in the Universe (Cappellaro et al. 1999). The 

. Gamma Ray Bursts (GRBs) are sudden and powerful flashes of classification can be further subdivide d in Types H-L and II-P, 

£ gamma-ray radiation occurring random at a rate of ~ 1 per day accordi to the sh of the H ht curve which can be linear (L) 

in the Universe. The duration of GRBs at MeV energies ranges Qr haye a , ateau (p) T ^ SNe haye jj lines in emis- 

, t0 ^ ° S6 f' g tS g character - sion. They can be distinguished from the classical Type II SNe 

ized by a duration > 2 seconds. b ^ lack of ab tion in their Balmer lines . | C hugai (1997) 



During the last years, it has been established that at least a introduced this designation to refl ect the fact that these SNe un 

large fraction of long-duration GRBs is directly connected with d si ificant intera ction with a "dense wind" produced by 

the death of massive stars. This scenario has got strong support ^ SN itor ior t0 eX pl os ion. Type la SNe are character- 

from observations of supernova (SN) features in the spectra of ized fe ^ lack of H Hnes and b a absorption observed 

several of GRB afterglows. Examples of the spectroscopic SN- . ,„.„ „„, s ^ , ^ ^_ _ „ ,.. £ „. TTrr 

„„„ .. • , j oh innou i™b nonmc /r i t at AA 6347, 637 1 A attributed to the P-Cyg profile of Si II. Type 

GRB connection include SN 1998bw/GRB 980425 (Galama et T . „. T , , Jb _ , 

i innoN cm onm JU pdd mmon /c, 1 t 1 onm u t u irj SNe are identmed by spectra with no evident Balmer lines, 
al. 1998), SN 2003dh/GRB 030329 (Stanek et al. 2003, Hiorth . , „ „. „ ,. J ^ , . TT T ,. tn — ^ , /.nfj i 

<- i cMinmi /r-DD rann ah i ■ <- i oaa/i\ a weak or absent Si II lines and strong He I lines. Bertola (1964) 

et al. 2003), SN 20031w/GRB 031203 (Malesam et al. 2004) and _^ , ~ „ , c 1 c o^t — i, . ^ t . T i... 

cm ™,: 'i™b (imoio . i ™,: \ i reported the hrst observation of this class of SNe, but the lb 

SN 2006ai/GRB 060218 (Masetti et al. 2006, Campana et al. , r . . . , , , . ^rr. . T 

?006 Modiay et al 9006 s ) designation was introduced later by |Elias et al. (1985J Type Ic 



SNe are characterized by weak or absent H and He lines and 
no evident Si II. They show Ca II H&K in absorption, the Ca II 



2006, Modjaz et al. 2006) 

Originally, SNe were classified in to two main types (I, II) on 

the basis of the presence (Type II) or absence (Type I) of H lines , r . . , . , „ ,.. . , 

. . r ,rr ■ , a. .i. \ t- i ■ . . near-IR triplet with a P-Cygm prohle, and O I in absorption. The 

in their spectra (Zwicky 1938). At the present time we distin- ,„ „ , * , , , — . ' TT — = 

. , r „ „ XT J , . j.., . . . Ic class was introduced by Wheeler & Harkness (1986) 

guish among hve SN types on the basis of the spectra obtained . i. ' — — — , 

at maximum light. Classical Type II SNe have prominent Balmer f Interestingly the SN features observed in the afterglows 

}V F of GRBs resemble those of Type Ib/c SNe (Delia Valle 2006; 

Send offprint requests to: E. Bissaldi, currently at Max-Planck Institut Woosley & Heger 2006, and references therein). In particular, 

fur Extraterrestriche Physik, GiessenbachstraBe 1, 85748, Garching, modeling of the SN lightcurves reveals that the SNe with GRB- 

Germany, e-mail: ebs@mpe.mpg.de connection have very large explosion energy and mass produc- 
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tion of 56 Ni compared to normal Type Ib/c SNe (Iwamoto et al. 
1998, Nakamura et al. 2001, Deng et al. 2005, Mazzali et al. 
2006b), except SN 2006aj which requires an explosion energy 
that is comparable to that of normal SNe Ib/c (Mazzali et al. 
2006a). These facts have motivated people to invent the term 
"hypernovae" for this special and much more powerful class of 
SNe (Iwamoto et al. 1998; see also Paczynski 1998). 

These phenomenological results greatly support the "col- 
lapsar" model (MacFayden & Woosley 1999, Zhang, Woosley 
& MacFayden 2003) in which a Wolf-Rayet progenitor under- 
goes core collapse, producing a rapidly rotating black hole sur- 
rounded by an accretion disk which injects energy into the sys- 
tem and thus acts as a "central engine". The energy extracted 
from this system supports a quasi-spherical Type Ib/c SN explo- 
sion and drives collimated jets through the stellar rotation axis 
which produce the prompt gamma-ray and afterglow emission 
(e.g., Zhang & Meszaros 2004). 

For what said before, it is clear that a comparison between 
the Type Ib/c SN rates in galaxies and the GRB rate is very im- 
portant in order to test if the SN-GRB connection is in agree- 
ment with the actual knowledge about SN progenitors, stellar 
and galactic evolution. 

Galaxies of different morphological type evolve in a differ- 
ent way with the main driver of their evolution being the star 
formation history. Recently, Calura & Matteucci (2006) calcu- 
lated the rates of Type II, lb and Ic SNe in galaxies of different 
morphological type (ellipticals, spirals and irregulars), by means 
of self consistent galaxy evolution models, reproducing both the 
chemical and spectro-photometric properties of galaxies. 

In this paper, we present a study of models which accurately 
reproduce the local Type Ib/c SN rate in galaxies, by following 
the approach of Calura & Matteucci (2006). We assume differ- 
ent progenitors for the Type Ib/c SNe such as single Wolf-Rayet 
(WR) stars and massive stars (12 - 2OM ) in close binary sys- 
tems, as well as several star formation histories typical of galax- 
ies of different morphological type. 

The assumed range for the binary progenitors of Type Ib/c 
SNe are based on suggestions by Baron (1992), Nomoto et al. 
(1996) and on the recent modelling of the proto-type Ic super- 
nova SN1994 1 (Sauer et al. 2006), which suggests that the upper 
mass limit for these binary stars should not be larger than 2OM . 
For what concerns the single WR stars as Type Ib/c progenitors, 
we need to assume a lower mass limit for the formation of WR 
stars but this limit in uncertain and strongly depends on the de- 
tails of stellar evolution models such as rotation and assumed 
mass loss, which in turn depends on the initial stellar chemical 
composition (see Meynet & Maeder (2005) and Yoon, Langer 
& Norman (2006) for a recent study on the single WR stars and 
GRBs). Therefore, here we will explore several cases where the 
lower mass for the formation of a single WR star varies from 25 
up 90 M Q . However, it is important to recall that most stars with 
M > 25 M probably produce black holes, either in a prompt col- 
lapse or by fallback; in neither case one would expect a normal 
Type Ib/c SN, since probably most of the mass in the Si shell, 
which is turned into Ni during the explosion, would fall into 
the black hole. So, it is still unproven that massive SNe forming 
black holes would appear as typical SNe Ib/c. If, instead, such 
massive stars are in interacting binary systems, their final core 
evolution is drastically changed and they are expected to pro- 
duce neutron stars rather than black holes (Brown et al. 2001) 
and hence possibly become SNe Ib/c. All these facts, seem to 
suggest that the binary channel for the formation of SNe Ib/c is 
favored over the channel of single very massive stars. 



We then derive the SN rates and compare them with the ob- 
served local GRB rates, thus obtaining an estimate of the local 
GRB/SN ratio. Furthermore, we move on cosmological scales 
by computing the cosmic Type Ib/c SN rates and, on the basis of 
that, we chose the best model for the GRB distribution as a func- 
tion of redshift and therefore we predict the expected number of 
GRBs at high redshift. 

The paper is organized as follows: in Section 2 we present 
the galactic models and the theoretical Type Ib/c SN rates for 
galaxies of different morphological type. In Section 3 we dis- 
cuss the local GRB rate and compare it with the theoretical pre- 
dictions. In Section 4 we describe how we compute the cosmic 
Type Ib/c SN rates and we extract a best model for the distribu- 
tion of GRB as a function of redshift. Finally in Section 5 some 
conclusions are drawn. 

1 .0.1 . Evolution Models for Galaxies 

In order to compute the SN Ib/c rates and compare them with the 
observed GRB rate we adopted self-consistent galaxy evolution- 
ary models already reproducing the majority of the properties 
of galaxies of different morphological type. In particular, these 
models are tuned to reproduce the present time SN rates and the 
chemical abundances in stars and gas. Moreover, the predicted 
star formation and chemical enrichment histories, when adopted 
in spectro-photometric models, give good agreement with pho- 
tometric properties of galaxies. 

In order to better understand the adopted galaxy models, a 
schematic outline of their basic equations and physical processes 
is provided (Calura & Matteucci 2003). 

Let G,(f) be the normalized fractional mass of gas within a 
galaxy in the form of the element i: 



Gi{f) 



M g (t)X,(t) 

M tnt 



(1) 



where M, , and M g (t) are the total galaxy mass and the mass of 
gas at time f, respectively. 



Xi(t) 



Gj(t) 
G(t) 



(2) 



represents the abundance in mass of the element i, the summa- 
tion over all elements in the gas mixture being equal to unity. 
Thus, the quantity 



G(t) 



M g (t) 

M tn , 



(3) 



is the total fractional mass of gas present in the galaxy at time t. 

The temporal evolution of G,(0 is described by the basic 
equation: 



G, = -MXM + Ri(f) + (Gi). nf - (G,) 



(4) 



Here, if/(t) is the star formation rate (SFR), namely the fractional 
amount of gas turning into stars per unit time. Rj(f) represents 
the returned fraction of matter in the form of an element ;' that 
the stars eject into the interstellar medium (ISM) through stellar 
winds and SN explosions. This term contains all the prescrip- 
tions regarding the stellar yields and the SN progenitor models. 
The two terms (G,) and (G,) account for the infalling ex- 

V /inf \ 'out 

ternal gas from the intergalactic medium and for the outflow, 
occurring by means of SN driven galactic winds, respectively. 

Here we are interested in computing the SN rates in the 
framework of these self-consistent galaxy models. The two key 
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ingredients for computing the SN rates are the SFR and the stel- 
lar initial mass function (IMF). The first is generally expressed 
as a function of time only, while the second, which describes the 
stellar mass distribution at birth, is likely to be universal (Kroupa 
2002) and not vary as a function of time (Chiappini et al. 2000). 

1.0.2. Star Formation Rates 

The main feature characterizing a particular morphological 
galactic type is represented by the prescription adopted for the 
star formation history, summarized in the SFR expression. 
In the case of elliptical galaxies, the SFR ifr(t) (in Gyr^ 1 ) has a 
simple form and is given by: 



m = vG(t) 



(5) 



The quantity v is the efficiency of star formation, namely the 
inverse of the typical time scale for star formation. 

In the case of ellipticals, v is assumed to drop to zero at the 
onset of a galactic wind, which develops as the thermal energy of 
the gas heated by supernova explosions exceeds the binding en- 
ergy of the gas (Arimoto & Yoshii 1987, Matteucci & Tornambe 
1987). This quantity is strongly influenced by assumptions con- 
cerning the presence and distribution of dark matter (Matteucci 
1992); for the model adopted here a diffuse (7^/7^=0. 1, where 
R e is the effective radius of the galaxy and Rj is the radius of the 
dark matter core) but massive {Mtj ar k/Mi um = 10) dark halo has 
been assumed. 

In the case of irregular galaxies, a continuous star formation 
rate is assumed as in eq. [3] but characterized by a lower effi- 
ciency than the one adopted for ellipticals. 

In the case of spiral galaxies, the SFR expression is slightly 
more complex, since it takes into account the feedback between 
SNe and ISM by means of a dependence upon the total surface 
mass density (see Chiosi 1980, Chiappini et al. 1997): 



<Kr,f) 



' cr(r,f) 


2(k-l) 


v(r, t C al) 


<r(r ,t) 




cr(r, t) 



k-\ 



G k (r, t) 



(6) 



where v is the SF efficiency, cr(r, f) is the total (gas + stars) 
surface mass density at a radius r and time t, and <x(r Q , t) 
is the total surface mass density in the solar region. For the 
gas density exponent k, a value of 1.5 has been assumed by 
Chiappini et al. (1997)| in order to ensure a good fit to the obser- 
vational constraints at the solar vicinity and in agreement with 
the estimates by |Kennicutt (1998)[ 

The three different star formation rates for elliptical, irregu- 
lar and spiral galaxy models as functions of time are shown in 
Figure Q] The SFR of the elliptical model (corresponding to a 
galaxy with luminous mass of 1O U M ) is characterized by very 
high SFR values (from 100 to 1000 M yr~') and by a starburst 
lasting ~ 0.2 Gyr, which traces the rapid collapse of a homoge- 
neous sphere of primordial gas where star formation is taking 
place at the same time as the collapse proceeds. Star formation 
halts as the energy of the ISM, heated by stellar winds and SN 
explosions, balances the binding energy of the gas. At this time a 
galactic wind occurs, sweeping away almost all the residual gas. 

The SFR of the irregular model is continuous and is char- 
acterized by a low SF efficiency. In this scenario, irregulars are 
assumed to assemble from merging of protogalactic small clouds 
of primordial chemical composition, until a mass of ~ 6 x 10 9 
M is accumulated, therefore producing stars at a lower rate than 
spiral and elliptical galaxies. 

The predicted SFR for spirals is characterized by two peaks 
due to two infall episodes present in the model of Chiappini et al. 



(1997). During the first episode, the halo forms and the gas shed 
by the halo rapidly gathers in the center, yielding the formation 
of the bulge. During the second episode, a slower infall of exter- 
nal gas forms the disk, with the gas accumulating faster in the 
inner than in the outer region ("inside-out" scenario; Matteucci 
& Francois 1989). The process of disk formation is much longer 
than the halo and bulge formation, with timescales varying from 
~ 2 Gyr in the inner disk to ~ 7 Gyr in the solar region and up 
to 15-20 Gyr in the outermost disk. The two infall episodes are 
clearly visible in Figure Q] The second one is dominated by the 
contribution of the disk. Here, the SFR for the MW represents 
an average between 2 and 20 kpc from the galactic center, and it 
is assumed as representative for an average spiral galaxy. 

1.0.3. Initial Mass Function 

The most widely used functional form for the IMF is an ex- 
tension of that proposed by Salpeter (1955) to the whole stellar 
mass range: 



SalpeteAM) - ^Salpeter M 



-(!+*) 



(7) 



where x = 1.35, tflsaipeter — 0.17 is the normalization constant 
derived from: 



I' 

J M, 



M< 



'Salpeter 



(M) dM = 1 



(8) 



The IMF in eq. (7) is by number with M ]ow = 0.1 M a and 
Mio W = 100 M Q . We are adopting this IMF for ellipticals and 
irregulars, where there is not a direct measure of the IMF. On 
the other hand, multi-slope expressions of the IMF give a better 
description of the luminosity function of the main sequence stars 
in the solar neighbourhood, as proposed by Tinsley (1980) and 
Scalo (1986, 1998). A simplified two-slope approximation to the 
actual Scalo (1986) IMF, that we adopt for the Milky Way and 
spirals, is given by: 



4>s 



ffl/o86 



where 



(M) = ttsc, 



nloSb 



M 



■a+x) 



(9) 



x = 1.35 &s caloS6 - 0.19 if M <2M G 
x = 1.7 yt Wo86 = 0.24 if M > 2M . 

Also for this IMF holds the normalization condition expressed in 
eq. (8) for the mass range 0. 1-100 M e . Finally, one can consider 
the so-called "top-heavy" IMF, which was suggested by obser- 
vational and theoretical arguments (e.g. Larson 1998 and refer- 
ences therein) and which assumes that the mass scale of the IMF 
was higher at earlier stages of the Universe, i.e. at high redshift: 



4>,- h (M) = 3Kt- h M 



•(l+x) 



(10) 



where x — 1.0 and J\ t -h 0.14. However, this kind of IMF 
does not produce satisfactory results for the solar neighbourhood 
and it is better to assume a constant IMF (Chiappini et al. 2000). 

1.1. SNRi b / c models 

Generally, Type Ib/c SNe occur in the vicinity of star forming re- 
gions and usual candidates for their progenitors are represented 
by massive WR stars (e.g. Hamuy 2003). However, the limiting 
mass for the formation of a WR star depends on the assumptions 
made in stellar models such as mass loss rate and rotation as re- 
cently shown by Meynet & Maeder (2005). Their results suggest 
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Fig. 1. Star formation rates, expressed in M yr , as functions of time calculated according to equations[5]and[6]for elliptical (upper 
curve), spiral (median curve) and irregular (lower curve) galaxy models. 



that the limiting mass for forming a WR star increases with de- 
creasing metallicity because of the growth of the mass loss rate 
with metallicity. For the sake of simplicity we will start here by 
assuming that all massive stars with initial mass M > 25 M are 
becoming WR stars and explode as type Ib/c SNe. Then we will 
test also larger values for this mass limit. 
The type Ib/c SN rate is then expressed as: 

SNR ModelI (f) = <f>(M) \p(t - Tm) dM 

J 25 

J- 100 
<f>(M)dM, (11) 
25 

where tm is the stellar lifetime. Hereafter, equation QT| will be 
referred to as Model I. 

However, some observational evidence seems to be against 
the WR progenitor hypothesis. In particular, Type Ib/c SNe can 
show broader light curves than expected assuming a massive star 
collapse origin. Other problems are represented by the apparent 
paucity of WR stars, unable to account for the rates observed 
in local galaxies (Muller et al. 1992), along with the fact that 
some events have been observed relatively far from active star 
forming regions (Filippenko & Sargent 1991). More promising 
candidates, which could account for these evidences, are repre- 
sented by less massive stars, with typical masses of 12-20 M 
(e.g. Baron 1992), in close binary systems. In this case, the loss 
of their envelope occurs by means of Roche-Lobe overflows. In 
this scenario (Model II), the SN rate can be expressed as: 

^20 

SNR ModelII (f) = F <p(M)4 f (t-T M )dM 
J 12 

,-20 

^Fi//(t) 4>(M)dM (12) 

J 12 

where F is a parameter and represents the fraction of massive 
binary systems of that particular kind giving rise to Type Ib/c 
SNe. This parameter is chosen to be equal to 0.15 (Calura & 



Matteucci, 2006). This choice is motivated by the facts that first, 
in any galaxy, half of the massive stars are in binary systems, and 
second, that the fraction of massive stars in close binary systems 
(namely those which can give rise to SNe Ib/c) is reasonably ~ 
30%, i.e. similar to the close binary frequency predicted for low 
mass systems (Jeffries & Maxted 2005). Therefore the estimated 
value for this parameter is given by: 

F a 0.5 ■ 0.3 0.15. (13) 

This is probably a lower limit for the fraction of binary 
systems giving rise to SNe Ib/c, as suggested in a paper by 
Podsiadlowski et al. (1992), where they concluded that 15-30% 
of all massive stars (M > 8M ), can produce SNe lb. Moreover, 
a recent estimate by Kobulnicky et al. (2006) seems to favor the 
higher estimate. However, our choice of the binary fraction is 
also tuned to reproduce the observed SN Ib/c rates in galaxies, as 
shown by Calura & Matteucci (2006), when also single WR stars 
are taken into account. Therefore, taking both scenarios into ac- 
count, a cumulative SN rate is obtained (Model III): 

SNR ModelIII (f) = SNR ModelI (f) + SNR Model II (f) 
/ /~20 plOO \ 

* ip(f) I F J <D(M) dM + J O(M) dMj . (14) 

In this case, the global fraction of Type Ib/c SN progeni- 
tors is similar to the above fractions of massive binaries by 
Podsiadlowski et al. (1992). The three SN rates (Models I, II 
and III) as functions of time in the local Universe, i.e. for a spi- 
ral galaxy like the Milky Way, are plotted in Figure [2] A Scalo 
(1986) IMF (eq. 9) is assumed. The two main infall episodes, 
which are characteristic for the halo and for the disk formation, 
respectively, are clearly visible. 

Observational rates measured by different authors are also 
indicated at present time (t ~ 13.47 Gyr). The SN rates are 
usually expressed in supernova units, SNus, (1 SNu = 1 SN per 
100yrperl0 I0 L BQ ; 1 SNuK = 1 SN per 100 yrper 10 10 L K0 ). 
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Fig. 2. SN rate for a spiral galaxy as a function of time calculated according to Model I (equationQj] dotted line); Model II (equation 
[121 dashed line); Model III (equation[T4l solid line). Data from Cappellaro et al. (1999) {open triangle), Mannucci et al. (2005) {open 
square), Delia Valle (2005) {open rhombus) and Podsiadlowski et al. (2004) {star). A Scalo 1986 IMF is assumed. 



The rate of Type Ib/c SN exploding in spiral galaxies, i.e. type 
Sbc/Sd, was measured by Cappellaro et al. (1999) obtaining: 



SNRib /c (f ~ 13.47 Gyr) = (0.14 + 0.07) SNu. 

For the same rate, |Mannucci e t al.(2005) calculated a value: 

SNRib /c (?~ 13.47 Gyr) = (0.067 SNuK. 

In order to compare these rates with the SN rate models, which 
are calculated in units of SNe yr -1 , one needs to multiply them 
by the correct galactic luminosities. For spiral galaxies, the 
Milky Way B-band luminosity L B = (1.8 ± 0.3) x 10 10 L B0 (van 
der Kruit 1986) and K-band luminosity L K 6.7 x 10 10 L KQ 
(Kent et al. 1991) are assumed. 

Furthermore, the validity of the choice for the adopted mas- 
sive binary star parameter, namely F = 0.15, has been checked, 
performing a fine tuning with different values 0.1 < F < 0.5 
(Figure |3), both for the rates calculated from Model III (cumu- 
lative model) and Model II (binary massive stars model). It is 
evident that F - 0.15 represents quite a reasonable choice and 
it is in good agreement with the observations, especially for the 
cumulative model (Mod III). This is consistent with the results 
by Calura & Matteucci (2006). 

Extending the analysis also to elliptical and irregular 
galaxies, the results for the cumulative SN rate model are shown 
in Figure [4] A Salpeter IMF has been adopted for modeling the 
SN rate in ellipticals and irregulars. F = 0.15 is assumed. The 
observed rates are indicated for comparison in the correct units. 
In the case of irregular galaxies, the conversion from SNus is 
calculated by applying the B-band luminosity of the main body 
of the Small Magellanic Cloud (Vangioni-Flam et al. 1980), i.e. 
V r B res - 7.8 x 10 s L B0 . 

Therefore, the rate measured in irregulars by 



Cappellaro etal. (1999) 



SNRib /c (f ~ 13 Gyr) = (0.22 + 0.16) SNu, 



is converted in to a value of: 

SNRib /c (f ~ 13Gyr) = (1.72 + 1.2) x 10~ 4 yr" 1 . 

^From Figure|4]it is evident that the SN rate model for an irreg- 
ular galaxy well reproduces this observed rate. 



1.2. The local GRB rate 

The rate of observed GRBs in a galaxy such as the Milky Way 
can be established from the BATSE monitoring as R b s - Sokolov 
(2001) calculated this value assuming that all or at least the 
main part of long GRBs were associated with SNe Ib/c. From 
the observations of known GRB host galaxies between 1997 and 
2000 (Bloom et al. 2002), he deduced that the absolute magni- 
tude MnostCai ranged from 22 to 28.5 mag. The search for direct 
GRB-SN associations in nearby galaxies represents a challenge, 
because the majority of the SNe related to GRBs are faint (22- 
26 mag) and in very distant galaxies with z > 0.4^1.5. Sokolov 
(2001) considered the number of galaxies brighter than 26 mag 
in one square degree of the sky, N ga ;< 26 mag ~ 2x 10 5 galaxy 
deg~ 2 (Casertano et al. 2000), and the number of GRBs observed 
by BATSE, ncRB ~ 0.01 deg~ 2 yr -1 . ^From these two quantities 
he calculated a rate of GRB events R Q b s ~ 5.0 x 10~ 8 galaxy^ 1 
yr -1 . 

Other local values calculated by different authors are reported in 
Table[[](/e/f table) in units of GRB yr" 1 . 



1.2.1. "True" GRB rates 

The so-called true GRB rate is given by: 

Rgrb = {ft 1 } Robs- (15) 

where the quantity f h l is the beaming factor (Sari et al. 1999), 
which accounts for the fact that a GRB does not light up the full 
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Fig. 3. SN rate for a spiral galaxy as a function of time calculated according to Model II (equation 12, Left Panel) and Model III 
(equation 14, Right Panel). Different values of the parameter F are considered: 0.1 (solid line); 0.2 [dotted line); 0.3 [short-dashed 
line); 0.4 (long-dashed line); 0.5 (dot-dashed-line); Data from Cappellaro et al. (1999) (open triangle), Mannucci et al. (2005) (open 
square), Delia Valle (2005) (open rhombus) and Podsiadlowski et al. (2004) (star)A Scalo 1986 IMF is assumed in both cases. 



o.i 



0.01 



0.001 



0.0001 



Ellipticals 



Spirals 



Irregulars 



0.001 



0.01 0.1 
Time [Gyr] 



10 



Fig. 4. Model III SN rate as a function of time for galaxies of different morphological type: Ellipticals (solid line), Spirals (dashed 
line) and Irregulars (dotted line). ,F=0.15 is assumed. Data from Cappellaro et al. (1999) (open triangle), Mannucci et al. (2005) 
(open square), Delia Valle (2005) (open rhombus) and Podsiadlowski et al. (2004) (star); experimental data for the local SN rate in 
irregular galaxies (open circle) are also shown. 



Table 1. 



Local GRB Rates R obs 


GRB Beaming Factor 


Author 


GRByr 1 


Author 


ft 1 


) 


Sokolov (2001) 

Zhang & Meszaros (2004) 

Firmani et al. (2004) 

Delia Valle (2005) 

Le & Dermer (2006) 


~ 5.0 x 10"* 
~ 4.0 x 10~ 7 
~ 2.0 x 10~ 7 
~ 3.8 x 10- 7 
~ 1.5 - 1.9 x 10" 5 


Schmidt (2001) 
Frail et al. (2001)| 


~ 100 
520 ± 85 
75 + 25 
~ 340 
13 < (f- 1 ) < 10 4 


Guetta et al. (2003) | 
fYonetoku et al. (2005)| 
[Soderberg et al. (2005) 
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Fig. 5. SN rate for a spiral galaxy as a function of time calculated according to Model I (equationQj] dotted line); Model II (equation 
[T2l dashed line); Model III (equation [T4l solid line). A Scalo 1986 IMF and F = 0.15 are assumed. Upper and lower limits for the 
observed local GRB rates (equationfTol) are over-plotted (upper and lower arrows). The open rhombus represents a "typical" value 
calculated from a local rate by Zhang & Meszaros (2004) and a beaming correction by uncorrected value by |Frail et al. (200Tj| (see 
equation [17). 



celestial sphere but rather a fraction. Since the possible correc- 
tions for the beaming factor cover a large range of values (see 
Table [1] right table), only plausible upper and lower limits can 
be estimated for the local GRB rate: 



R 



GRB, min 



R 



GRB, max 



(■4, min) ' R°bs, 

3 x lO^yr" 1 

(/fc.max) ' R 

4 x 10~ 4 yr 



obs, max 
-1 



(16) 



where only values calculated after 2004 were considered. 

A "typical" value that is often considered in literature is 
given by the local rate calculated adopting the uncorrected 
value by Zhang & Meszaros (2004) and the beaming factor by 
(2001)| namely: 



|Frail et ; 



(17) 



R GRB - (520 ± 85) • (4.0 x 10~ 7 ) yr 1 

~ (2.08 + 0.34) ■ 10~ 4 yr" 1 . 

Figure [5] shows the distributions of the three SN rate Models 
(from equations [TTl[T2l and[T4l) for a spiral galaxy comparedwith 
the observed local GRB rates. Again, a Scalo (1986) IMF and 
F = 0.15 are assumed for Model I— III. The arrows mark the 
lower and upper limits for the GRB rates as calculated in equa- 
tions [16] while the open rhombus represents the "typical" GRB 
rate calculated in equation[T7l 

As one can see from Figure [5] the SN rate Models poorly 
reproduce the observed local GRB rates. However, uncertainties 
regarding Type Ib/c SN progenitors masses allow us to perform 
a tuning of Model I (massive single star model, see equation 
ITTb . In particular, we considered different mass ranges for the 
progenitors of Type Ib/c SNe, namely: Ami = 25-100 M ; Am2 
= 40-100 M Q ; A M3 = 60-100 M G ; A M4 = 80-100 M ; A M5 = 90- 
100 M . The behaviours of the SN rate calculated with Model 



I and Model III in spirals are then shown in Figure [6] (left and 
right panel, respectively). This is allowed by the fact that the 
limiting mass for the formation of a WR star is still uncertain 
and depends on the stellar mass loss, which in turn depends on 
the stellar metallicity. 

Taking a closer look to Model I, it is evident that the range of 
observed GRB rates, in particular the one between the "typical" 
value (open circle) and the lower limit (bottom arrow), is better 
reproduced if M;„/ > 40 M Q is assumed. This is consistent with 
the observations of the four GRBs with a clear SN association 
(see Introduction), which predict a progenitor mass as large as 
40 M Q in three out of four cases. 

A last interesting distribution to study is the SN rate for 
Model III in galaxies of different morphological types together 
with the observed local GRB rates (Figure[7]i. Again, a Salpeter 
IMF is adopted for ellipticals and irregulars. Here the integration 
is performed over the mass range Am2 = 40-100 M for single 
massive stars giving rise to WRs. The parameter F = 0.15, as 
previously discussed. It is evident that in all cases the SN rate is 
in good agreement with the observed local GRB rates. It is worth 
noting that in the case of irregular galaxies, the predicted SN rate 
falls exactly in the range between the "typical" local GRB rate 
and its lower limit, which is consistent with the latest observa- 
tions of GRB-SN hosts (e.g., Conselice et al. 2005; Savaglio et 
al. 2006) and with a recent work by Fruchter et al. (2006) , who 
showed that the host galaxies of GRBs appear to be significantly 
faint and most of them are irregulars. 



2. Cosmic SN Ib/c and GRB rates and GRB/SNe 
ratios 

The comoving rate density is a measure of the number of events 
occurring per unit comoving volume and time which provides 
a census of the number of objects formed at a given redshift 
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Fig. 6. SN rate calculated for Model I (Left Panel) and Model III (Right Panel) in a spiral galaxy for five different mass ranges of 
existence of Type Ib/c SNe progenitors (single WR stars): Ami = 25-100 M Q (solid line); Am2 = 40-100 M (dotted line); Am3 = 
60-100 M (short-dashed line); A M 4 = 80-100 M (log-dashed line); A M 5 = 90-100 M Q (dot-dashed line). A Scalo (1986) IMF 
and a value of F = 0.15 are assumed. Upper and lower limits for the observed local GRB rates (equation[T6l> are over-plotted (upper 
and lower arrows). The open rhombus represents a "typical" value calculated from a local rate by Zhang & Meszaros (2004) and a 
beaming correction by uncorrected value by Frail et al. (2001) (see equationfTTl). 
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Fig. 7. Upper and lower limits for the observed local GRB rates (upper and lower arrows) compared with Model III SN rate 
model for galaxies of different morphological type: ellipticals (solid line), irregulars (dotted line) and spirals (dashed line). The 
open rhombus represents a "typical" value calculated from a local rate by Zhang & Meszaros (2004) and a beaming correction by 
uncorrected value by |Frail et al. (20 01) (see equation[T7b. 



Table 2. Values of the coefficients for i//mdp, (0> <Amdp 2 (0 an d iAstimW- 



SFR | 


a 


b 


c 


d e 


d 


P*MDP! (0 


0.049 


5 


0.64 


0.2 




P«MDP 2 W 


0.336 


5 


0.64 


0.0074 0.0197 


1.6 


P»STMw 


0.182 


1.260 


1.865 


0.071 
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and which can help understanding the object/structure formation 
in its various stages of evolution. In the following sections this 
quantity will be estimated and discussed for both Type Ib/c SNe 
and GRBs. 



2.1. The cosmic SN Ib/c rate 

In order to reproduce cosmic SN Ib/c rate densities, the funda- 
mental ingredients are a cosmic SFR density and a suitable IMF. 
There are different models which can describe the global SFR 
per unit comoving volume. In the following, four models will be 
reviewed and analysed. 



2.1 .1 . Cosmic SFR density models 



One model was computed by |Calura & Matteucci (2003)j (here- 
after CM03). They calculated the cosmic SFR density p* as a 
function of redshift (Madau's plot) in the case of a ACDM cos- 
mology and galaxy formation at Zf = 5. Since the cosmic SFR 
density is not a directly observable quantity, CM03 evaluated it 
by calculating the luminosity density (LD), p A , at certain wave- 
lengths and by adopting a universal IMF (the Salpeter one). 

The total LD in a given band is the integrated light radiated 
per unit volume from the entire galaxy population and it needs 
the luminosity function (LF) to be determined. The LF repre- 
sents the distribution of absolute magnitudes for galaxies of any 
specified Hubble type and it is often parametrized according to 
the form defined by Schechter (1976): 



(D(L) 



dL 



4>* 



exp 



L \ dL 



(18) 



where is a normalization constant related to the number of 
luminous galaxies per unit volume, L* is a characteristic lumi- 
nosity and y is the slope of the luminosity function. The LD 
stems from the integral over all magnitudes of the observed LF: 



PL 



I 



<D(L/L*)L/L* d(L/L±) . 



(19) 



The determination of the SFR density is related to the mea- 
sures of star formation in galaxies. CM03 reconstructed the 
history of cosmic star formation in the Universe by means 
of detailed chemical evolution models for galaxies of differ- 
ent morphological types, as described before. They normalized 
the galaxy populations to the B-band LF observed in the local 
Universe and studied the redshift evolution of the LD in various 
bands (U, B, I and K), calculating galaxy colors and evolution- 
ary corrections by means of a detailed synthetic stellar popula- 
tion model (Jimenez et al. 1998). The cosmic SFR density was 
computed according to: 



P*(z) = Y/>bj(z) (-£-) <Ai 



(z) 



(20) 



where p B j represents the B-band LD, yj;) B , * s tne ^" rjanc l mass- 
to-light ratio and i/r, represents the star formation rate for the 
galaxies of the 2th morphological type. 

An analytical fit which reproduces the total cosmic SFR den- 
sity as a function of time (CM03), is given by: 

y(t) = (0.337 -9.37 -t+ 11. 19 • f 2 ) 

for t< 0.328 Gyr; 



y(f) = (- 



1.57 + 0.11 -t + 0.05 



for 0.328 < t < 2.150 Gyr; 

y(t) = (-0.836 - 0.13 ■ t + 0.005 • ? 2 ) 

for t > 2.150 Gyr, where y(t) = p\. 

The SFR density for this model is then calculated as: 

Pcmo3(0 = lO v( ' ) M yr- I Mpc- 3 . 



(21) 



(22) 



The CM03 model predicts a peak at the redshift of galaxy forma- 
tion due to starbursts in spheroids. This peak is clearly visible at 
z ~ 5 (Figure 8 of CM03) and it is followed by a flat behaviour 
between z ~ 4.2 and z ~ 3 due to star formation in spiral galax- 
ies. The maximum star formation in spirals causes a smaller peak 
at z = 2; these galaxies are also responsible for the decline of 
the SFR density between z = 2 and z — 0, in agreement with 
observational data. 

Other two cosmic SFR density models considered for 
this analysis were calculated following the work done by 
Madau, Pozzet ti~& Dickinson (1998")] in which the emission his- 
tory of field galaxies is modeled at UV, optical and near-IR 
wavelenghts by tracing the evolution with cosmic time of their 
luminosity density, 



Pv(z) 



Jo 



L v y v (L v ,z)dL v 



(23) 



where y (L v , z) is the best-fitting Schechter luminosity function in 
each redshift bin. The integrated light radiated per unit volume 
from the entire galaxy population is an average over cosmic time 
of the stochastic, possibly short-lived star formation episodes of 
individual galaxies, and follows a relatively simple dependence 
on redshift. Madau, Pozzetti & Dic kinson (1998)] used a stellar 
evolution model defined by a time-dependent star formation rate 
per unit volume, i/r(f), a universal IMF and some amount of red- 
dening. In such a system, the luminosity density at time t is given 
by the convolution integral: 



Pv(t) = p e 



f 

Jo 



l v {t')i]/(t - t')dt' , 



(24) 



where l v (f) is the specific luminosity radiated per unit initial 
mass by a generation of stars with age l v (t'), p esc is a time- 
independent term equal to the fraction of emitted photons which 
are not absorbed by dust and the cosmic SFR density is derived 
from the observed UV luminosity density. 

Based on these considerations, the cosmic SFR 
model is given by an analytical fit developed by 
|Madau, Delia Valle & Panagia (1998)| (hereafterMDP98): 

P«mdp, (0 = a [ t b 9 e-'^ + d(l- e-^ c )] (25) 

expressed in M & yr~ l Mpc~ 3 and where is the Hubble time in 
Gyr, and the values for the coefficients are given in Table [2] 

This cosmic SFR was adopted by MDP98 in order to predict 
the evolution of the observed comoving luminosity density in 
the 'hierarchical clustering' scenario, in which elliptical galax- 
ies form continuously from the merger of disk-bulge systems or 
other ellipticals and most galaxies never experience star forma- 
tion rates in excess of a few solar masses per year. 

A second analytical fit by MDP98 consists in a large star for- 
mation density at high redshifts and was designed to mimic the 
'monolithic collapse' scenario as in the CM03 model, in which 
spheroidal systems form early and rapidly, experiencing a bright 
starburst phase at high-z: 



P*MDP,( f ) 



ae-' 9/f + d(l 



(26) 
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Fig. 8. Cosmic SFR densities as a function of redshift calculated for different epochs of galaxy formation, namely Zf = 3 (left panel) 
and Zf = 5 (right panel); Model CM03 (equation l22t: solid line; Model MDPi (equation |25li: short-dashed line; Model MDP2 
(equation|26|>: long-dashed line; Model ST03 (equation|27|>: dotted line. 
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Fig. 9. Same as Figure [8] calculated for different epochs of galaxy formation, namely zj — 8 (left panel) and zj = 10 (right panel). 
Models CM03, MDPj, MDP 2 , ST04 are indicated as in Figureg] 



expressed in units of M Q yr 1 Mpc 3 , where the values of the 
coefficients are given in Table [2] 

The last model which will be considered in the following 
analysis, was computed by Stolger et al. (2004, hereafter ST04), 
assuming a modified version of the parametric form of 41(f) as 
suggested by Madau et al. (1998). An analytical expression for 
this model is given by: 

PWO = a(t b e-" c + de«->°V c ) (27) 

always in units of M yr _1 Mpc -3 , where t is in Gyr. 
By fitting the measurements of SFR from several surveys 
|Giavalisco et al. (2004)| ST determined the coefficients of the 
function as summarized in Table[2] This model takes the correc- 



tions for extinction into account. Here t is the age of the Universe 
and ?o = 13.47 Gyr corresponds to z = 0. 

All the already described models are plotted together as 
functions of redshift in Figures [8] and [9] assuming four different 
epochs of galaxy formation Zf, namely 3, 5, 8 and 10. The differ- 
ent behaviours of the cosmic SFR following the 'monolithic col- 
lapse' scenario (Models CM03 and MDP2) or the 'hierarchical 
clustering' scenario (Model MDPi and ST04) are clearly visible 
toward higher redshifts. 

2.1.2. Cosmic SNib/c models 

The strategy for estimating the cosmic SN rate is the same as 
the one adopted for the local case in Section [TTTl Equations QT| 
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Fig. 10. Cosmic SN rate for Model III based on the cosmic SFR by CM03 and different IMFs: Salpeter (1955) (equation[7] solid 
line); "top-heavy" (equationfTUl dotted line); Tinsley (1980) (short-dashed line); Scalo (1986) (equation[9] long-dashed line); Scalo 
(1998) (dot-dashed line). Redshift of galaxy formation: z/ = 10. 



fT2l and [T4l for Model I— III, respectively, can be calculated con- 
sidering different cosmic SFRs, IMFs and various mass ranges. 
^From now on, the analysis will focus on the cosmic SN rate 
calculated for Model III, i.e. taking both massive single stars 
and massive binary stars into account, with a SFR by CM03. 
This is in fact the Model which gives the best fit to the assumed 
Type Ib/c SN rates in galaxies (Section 1.1). Moreover, the SFR 
will be computed assuming galaxy formation taking place at 
z f = 10. 

This choice of zj is motivated by the observation of the 
most distant quasars (Staguhn et al. 2005), suggesting mas- 
sive starbursts up to z ~ 5, together with an interesting result 
by Mobasher et al. (2005). Searching for high-redshift /-band 
"dropout" galaxies in the portion of the Great Observatories 
Origins Deep Survey (GOODS) southern field that is covered 
by extremely deep imaging from the Hubble Ultradeep Field 
(HUDF), they found evidence for a massive galaxy, HUDF- 
JD2, whose spectral energy distribution shows distinctive fea- 
tures which are consistent with those of a galaxy at z ~ 6.5, 
observed several hundred million years after a powerful burst of 
star formation. The best-fitting models adopted by Mobasher et 
al. placed the formation of the bulk of the stars at z > 9. For com- 
pleteness, they also reported alternative models of dusty galaxies 
at z ~ 2.5. However, these models provided significantly poorer 
fits to the photometric data. Mobasher et al. suggested that, if 
the high-redshift interpretation is correct, HUDF-JD2 is an ex- 
ample of a galaxy that formed by a process strongly resembling 
traditional models of 'monolithic collapse'. 

An interesting test is to model the SN rate for different IMFs. 
Figure [TO] shows its behaviour as a function of the redshift, for 
various choices of the IMF, namely a Salpeter IMF; a "top- 
heavy" IMF; a Tinsley (1980) IMF; a Scalo (1986) IMF, and 
a Scalo (1998) IMF. The Salpeter IMF, which lies at an interme- 
diate position among the different models, is chosen as the best 
IMF for performing further analysis. 



2.2. The cosmic GRB rate 

The effort to find a redshift estimate for bursts without optical 
afterglows represents an open research branch in the GRB cos- 
mology field. The redshift distribution, together with the LF, can 
provide important insights not only into the physics of the in- 
dividual objects themselves, but also into the evolution of the 
matter in the Universe. 

Often, when doing large statistical studies of a particular 
class of objects, the LF and redshift distribution are assumed to 
be independent quantities; that is, the luminosity function of the 
sources are assumed to be the same for all redshifts. This makes 
the analysis easier when one has limited information (e.g. signif- 
icant data selection effects) obscuring a direct interpretation of 
the measured distributions of luminosity and redshift. However, 
it has been shown that this assumption is not valid for many 
astrophysical objects, i.e. quasars (Boyle 1993), as well as for 
GRBs. 

By studying the prompt and especially the afterglow emis- 
sion in great detail, the GRB community realized from the very 
first bursts with measured redshifts that the GRB LF was not 
very narrow and in fact it exhibited a rather large dispersion, 
hence preventing the flux of a GRB to be used as a standard can- 
dle from which to infer a redshift for the source. 

In the last years, various types of "standard-candle relation- 
ships" for GRBs have been discussed, from which a redshift 
can be inferred based on common GRB observables. For exam- 
ple, Lloyd-Ronning, Fryer and Ramirez-Ruiz (2002, hereafter 
LFR02) studied GRB formation rates as the result of a redshift 
distribution estimated by Fenimore & Ramirez-Ruiz (2000), by 
means of a "Variability-Luminosity" correlation. Using a r sta- 
tistical method (Lynden-Bell 1971, Efron & Petrosian 1992), 
LFR02 obtained a GRB comoving rate density as: 

p™(z) =Po| (1+zr t for z > 2? (28) 
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where po is the normalization constant. This distribution was 2.3. Cosmic SA/u,/ c vs. GRB rates 



computed in units of Mpc~ 3 yr~' with an arbitrary normaliza- 
tion of the curve to ~ 1 at ( 1 + z) = 1 . 

Another attempt in the GRB redshift estimate was done 
by Yonetoku et al. (2004, hereafter Y04), adopting a "Peak 
Luminosity-Peak Energy" correlation. Introducing the same sta- 
tistical method as LFR02, they produced the following GRB rate 
density: 



Pyo4(z) = po 



(1 + z)6 .0±1.4 for 
(1 + z )0-6±0.2 for 



Z < 1 

z > 1. 



(29) 



In a recent paper by Matsubayashi et al. (2005, hereafter 
M05), the apparent cosmic GRB rates by LFR02 and Y04 have 
been studied and used to derive the true absolute GRB forma- 
tion rate. The peculiarity of this work lies in the fact that M05 
took into account the geometrical correction of the jet opening 
angle, 9j, and the jet-luminosity evolution, finally presenting an 
analytical formula to calculate the true GRB rate. 

M05 adopted the cosmic GRB formation rate proposed by 
Y04, performing the normalization of Pyoa(z) at z — 1, and ob- 
tained a true comoving GRB rate density given by: 



Pmo5(z) = Proi(z)(f b ! ) 



(30) 



They derived the distribution of the beaming factor (f b l \ in a 
very detailed way, obtaining: 



{fb% = 200 



l+z 



it-p 



(31) 



Here a and /3 are the slope indexes of the GRB luminosity evo- 
lution A(z) (see LFR02 and Y04 for a detailed description), 



A(z) oc (1 +zf, 



(32) 



and the jet-corrected luminosity evolution Lj(z) (Lamb et al. 
2005), 



Lj(z) oc (l+z) 



(33) 



respectively. M05 imposed 1 ) = 200 at z = 1 . In this 
case, the mean jet opening half-angle, defined by (o/j = 

[(■4 ' ) /2] ' > is equal to the typical value of 0. 1 rad at z = 1 . 

The value of a has been studied by many authors. LFR02 
found a =* 1.4 + 0.5, while Y04 calculated a 2.60+°^, 



which is quite similar to that of quasars. Firmani et al. (2004) 
found evidence supporting an evolving LF where the luminosity 
scales as a 1.0 + 0.2. M05 adopted the value by Y04, i.e. 
a = 2.6 in his work. 

Studying possible values for /3, M05 considered three cases: 
a) j8 = a - no jet-opening angle evolution; b) /3 = a/2 - in- 
termediate; c) /3 — - nojet-luminosity evolution. M05 rejected 
cases (a) and (c) (see the original work for an exhaustive descrip- 
tion), concluding that intermediate values of ft in the neighbour- 
hood of /3 - a/2, for which both Lj and (f^f depend on z, may 
be preferable. Applying a generalized r statistical method, M05 
calculated /3mo5 = 2.05^] ^. Recently, the existence of open- 
ing angle evolution of the form oc (1 + z) r with y = — 0.45*qjS 
was suggested by |Yonetoku et al. (2005)| Combining this with 
the result of a by Y04, a value of /3yo4 = a + 2y ^ 1.7+° ^ is 
obtained. 



Summing up all the results gathered up to this point, a compar- 
ison between the cosmic SN rate models and the cosmic GRB 
rate models is the last important step to make in order to obtain a 
complete overview. As it was explained in the previous Sections, 
the cosmic SN rate can be calculated in different scenarios and 
for different values of the parameters. Here, the cosmic SFR 
by CM03 is considered ('monolithic collapse scenario'), calcu- 
lated for a redshift of galaxy formation Zf = 10, and in units 
of M Q Gpc 3 yr _1 . Two different IMFs are assumed, namely 
the Salpeter IMF and the "top-heavy" IMF (see Section [TTX3l >. 
Model III is computed assuming the usual mass range Ami = 25- 
100 M and a massive binary parameter F = 0.15, as previously 
discussed. 

These SN rate models are compared to the cosmic GRB rate 
density models by LFR02, Y04 and M05. Models by LFR02 and 
Y04 need to be normalized to the observed local GRB rate, this 
time expressed in units of Gpc 3 yr . A typical value which is 
found in the literature is po = p(z = 0) = 0.5 Gpc 3 yr _I 
(Schmidt 2001). Taking different possible beaming corrections 
into account (see Table [U, a range for the true local GRB rate is 
obtained: 

p , mi „ ~ 38 ± 13Gpc- 3 yr _1 (34) 
po.max ~ 260 + 40Gpc- 3 yr- 1 . 

Figure [TT] shows the behaviours of the cosmic SN rate (up- 
per solid and dashed curves, for models by CM03 and ST04, 
respectively) and the estimated comoving GRB rate densities 
(shaded areas for models by LFR02 and Y04, and lower solid 
line for the model by M05). The SN rate in the left and right 
panels are calculated assuming a Salpeter and a "top-heavy" 
IMF, respectively. The GRB models by LFR02 (p LF m 2 (z)) and 
by Y04 (pyqa(z)) are normalized in the interval between po,min 
and po.max, corresponding to the (dark grey and light grey ar- 
eas, respectively. The GRB model by M05 (pmos(z), lower solid- 
line) is calculated according to the original parameters a = 2.6, 
/3 — 1.7 and p(0) 0.01 Gpc 3 yr -1 . The predicted cosmic 
SN rates are by far larger with respect to the estimated cosmic 
GRB rate models. However, this fact is not surprising. First, lo- 
cal GRB rates suffer from great uncertainties, mainly due to the 
beaming corrections. Moreover, the models applied for the cos- 
mic GRB rate densities rely upon luminosity indicators and con- 
sequent redshift estimates, which are still far from being con- 
firmed as "true" standard-candle relationships. It is worth not- 
ing that Figure 1 1 shows that the cosmic SN Ib/c rate obtained 
by means of the hierarchical scenario (Model ST04) predicts a 
much lower GRB rate at high redshift than the one predicted by 
the monolithic scenario (Model CM03). 

Looking again at FigureQT| one can see by eye that this ratio 
GRB/SNe is almost ~ 10~ 2 4- 10~ 3 , depending on the choice of 
the different parameters, which play an important role in the cal- 
culation of the SN rate. Table[3]summarizes in four sub-tables all 
the different analyses which have been performed in this paper, 
giving an estimate of the local ratio between GRB and Type Ib/c 
SNe according to the various applied models and parameters. 

In the tables, two ratios are presented, namely: 

Pi), max . m 

= 7777^ (Z = 0) (35) 



Kb = 



SNRn,/ c 

p0,min 



(z = 0), 



SNRi b/c 

i.e. the ratio are calculated assuming the largest and the smallest 
(fb 1 )■ This is to underline the great uncertainty of the local GRB 
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1 + z 1 + z 

Fig. 11. Absolute cosmic SN rate for Model III calculated for the SFR by CM03 (upper solid curve) and by ST04 (dashed curve), 
assuming first a Salpeter IMF (equation|7] left panel) and then a "top-heavy" IMF (equation[10] right panel). The plots show also 
three absolute comoving GRB rate densities, calculated according to the models by LFR02 (plfrqi(z), equationEFt. by Y04 (jdy(h(z), 
equation l29il and M05 (pmos(z), equation[30] lower solid line). Models by LFR02 and Y04 are normalized in an interval between 
two local GRB rates, po.min (bottom lines) and po.max (upper lines, see equation [34b. covering two shaded areas of possible values 
(light grey, superior and dark grey, inferior areas, respectively). 



rate due to the beaming correction. Still, the results are more or 
less consistent with previous works (Podsiadlowski et al. 2004, 
Delia Valle 2005, Le & Dermer 2006). 

While the scatter of possible values for the GRB-SN ratio 
is mainly due to uncertainties affecting the observed local GRB 
rate and the beaming factor, still the ratio is intrinsically small. 
This is motivated by the fact that only a very tiny fraction of 
all massive stars are capable of producing GRBs. Type Ib/c SN 
represent "obvious" candidates because they have lost the H en- 
velope when the collapse of the core occurs, therefore allowing 
the ultra-relativistic jets to escape from the massive progenitor. 
However this fact does not seem to be sufficient. According to 
the SN and GRB rates and to the estimates, only 1% - 

0.1% of SNe/i/ e are able to produce GRBs. This fact implies 
that GRB progenitors must have some other special character- 
istic other than being just massive stars. Recently, several au- 
thors have discussed a number of promising mechanisms at play 
in the GRB phenomenon, such as the stellar rotation (Woosley 
& Heger 2006), binarity (Podsiadlowski et al. 2004, Mirabel 
2004), asymmetric explosions (Maeda et al. 2006) and metal- 
licity (Fruchter et al. 2006, Langer & Norman 2006). 

3. Conclusions 

In this paper we have compared the local GRB rate with local 
Type Ib/c SN rates in galaxies of all morphological type: ellipti- 
cals, spirals and irregulars. To do that we have adopted for each 
galaxy type a specific history of star formation already tested in 
self-consistent models of galactic chemical evolution, which re- 
produce the main chemical and photometric properties of galax- 
ies and in particular the observed SN rates. We have assumed 
that Type Ib/c SNe, which are associated to long-duration, soft 
GBRs, originate either from massive single stars (M > 25M , 



WR stars) or from massive stars in the range 12-20 M , which 
are members of close binary systems. We made several tests by 
varying the IMF and the mass range for the existence of WR 
stars. This is justified by the still existing uncertainties in the 
progenitors of WR stars. 

Moreover, we computed the cosmic star formation rate, 
namely the star formation rate in a comoving unitary volume 
of the Universe, by taking into account all the galactic types and 
adopting a technique developed by CM03. Then the cosmic star 
formation rate has been used to compute the cosmic Type Ib/c 
SN rate as a function of redshift. In this particular formulation, 
which is the only one distinguishing among different galactic 
morphological types, we assumed that elliptical galaxies, espe- 
cially the most massive ones, form quickly and at high redshift, 
whereas the spirals and irregulars evolve much more slowly. The 
main driver of galaxy evolution is assumed to be the star forma- 
tion rate, being very fast in ellipticals and spheroids in general, 
and then decreasing going to spirals and irregulars. This sce- 
nario is different from the classic hierarchical scenario of galaxy 
formation, where ellipticals are assumed to have assembled by 
mergers and preferentially at low redshift. Because of this, we 
have tested also different cosmic star formation rates (and Type 
Ib/c SN rates) including the hierarchical scenario. We have then 
derived the GRB rate as a function of redshift and compared it 
with all the cosmic Type Ib/c SN rates considered. 

Our main conclusions can be summarized as follows: 

- The best agreement between the observed and predicted lo- 
cal Type Ib/c SN rate in galaxies of different morphological 
type is obtained by assuming that Type Ib/c SNe are both sin- 
gle massive WR stars and massive stars in binary systems, 
in agreement with previous results by Calura & Matteucci 
(2006). Therefore, we conclude that a model including both 
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Table 3. GRB/SN ratios R4 and Rg (equations [35j calculated according to the three models previously discussed (see equations 11, 
12 and 14). TABLE A displays Model I— III with the usual choice of the parameters, i.e. Am, = 25-100 M , F — 0.15 and a Salpeter 
IMF. Model III is our best model. TABLE B focuses on the SN rate calculated for the massive single star model, Model I, assuming 
4 different mass ranges (Am 2 -&m 5 ) and a Salpeter IMF, while TABLE C highlights the SN rate calculated for the massive binary star 
model, Model II, calculated for 5 different values of the binary parameter F and for a Salpeter IMF. Finally, TABLE D displays the 
cumulative model, Model III, calculated assuming the usual parameters, Am, and F — 0. 15, but different assumptions for the IMF. 



TABLE A 





Ra 


Rb 


Model I ,a 


~ 1.2 X 10~ 2 


~ 9.7 x 10~ 4 


Modelll (b c) 


~ 4.9 x 10~ 2 


-4.1 x 10~ 3 


Model III (a b c) 


~ 9.4 x 10" 3 


~ 7.9 x 10" 4 


TABLE B 


Model I c| 


Ra 


Rb 


A M , = 40-100 M Q 


~ 2.6 x 10~ 2 


~ 2.2 x 10~ 3 


A M ~ = 60-100 M 


~ 6.5 x 10- 2 


~ 5.4 x 10- 3 


A M4 = 80-100 M e 


~ 1.8 x 10"' 


~ 1.5 x 10~ 2 


A M5 = 90-100 M 


~ 4.2 x 10" 1 


~ 3.5 x 10" 2 


TABLE C 


Modelll lc) 


Ra 


Rb 


F = 0.1 


~ 7.4 x 10~ 2 


-6.1 x 10~ 3 


F = 0.2 


~ 3.7 x 10~ 2 


~ 3.1 x 10~ 3 


F = 0.3 


~ 2.5 x 10~ 2 


~ 2.0 x 10~ 3 


F = 0.4 


~ 1.8 x 10~ 2 


~ 1.5 x 10~ 3 


F = 0.5 


~ 1.5 x 10" 2 


~ 1.2 x 10" 3 


TABLE D 


Model III (a b) 


Ra 


Rb 


"top-heavy" IMF 


~ 3.2 x 10~ 3 


~ 2.7 x 10~ 4 


Tinsley IMF 


~ 1.6 x 10~ 2 


~ 1.3 x 10~ 3 


Scalo 1986 IMF 


~ 2.2 x 10~ 2 


~ 1.9 x 10~ 3 


Scalo 1998 IMF 


~ 8.6 x 10" 3 


~ 7.2 x 10" 4 



(fl| A M , = 25-100 M, 
m F = 0.15 
w Salpeter IMF 



progenitors should be preferred to compare with the GRB 
rate. 

- When comparing the local GRB rate with the predicted lo- 
cal Type Ib/c SN rate, the best agreement is reached under 
the assumptions that WR stars originate in the mass range 
40- 1OOM , and the predictions for a typical irregular galaxy 
are those which best fit the local GRB rate. This seems to 
confirm the recent finding by Fruchter et al. (2006) indicating 
that the hosts of GRBs are mainly irregular galaxies. Since ir- 
regular galaxies are normally metal poor objects, this would 
explain the need of having WR stars only for M > 4OM as 
due to the generally less efficient mass loss for metal poor 
stars (see Maeder (1992); Meynet & Mader (2005); see also 
the most recent papers by Yoon & Langer (2005) and Yoon et 
al. (2006) for the connection between WR stars and GRBs.) 

- The predicted cosmic SN rate computed with the cosmic star 
formation rate of CM03 follows the derived GRB rate as a 
function of redshift and predicts more GRBs at high redshift 
than hierarchical models. In particular, the predicted rate is 
larger that the cosmic GRB rate, for the whole lifetime of 
the Universe. This is not surprising since very likely only a 
small fraction of Type Ib/c SNe will give rise to GRBs. 



- The ratio between cosmic Type Ib/c SN rate and GRB rate, 
as predicted by our best model, is found to be in the range 
of si 10~ 2 - 10~ 3 for a large varieties of cases (different 
IMFs, different assumptions on Type Ib/c SN progenitors). 
We conclude that only 0.1 - 1% of SNe Ib/c are able to 
produce GRBs. This means that other special characteristics 
must play a role in creating GRBs: stellar rotation, asymmet- 
ric explosions and metallicity are all possible factors. 

- Because of the high star formation rate in spheroids at high 
redshift assumed in our best models, we predict a higher 
number of GRBs at high redshift than hierarchical models. 

Acknowledgements. E. B. wish to thank Massimo Delia Valle for interesting 
discussions and many stimulating comments and an anonymous referee for his/ 
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